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I. INTRODUCTION 

The quantum nature of gravity is not yet cleared in spite of endeavor of many researchers. 
An old idea on this issue is that gravity emerges as quantum effects of matter fields [1]. 
Originally, in such an induced gravity scenario, the Newton constant is naturally obtained 
from the one-loop calculation with a cutoff of the Planck scale. In this case, the induced 
cosmological constant becomes a huge amount if no special choice of the matter-field content 
is considered. 

We will consider a calculable model for induced gravity in the present paper. For this 
purpose, we first fix the choice of matter species to cancel the UV divergences. Next we 
should consider the mass spectra of the fields, which affect the finite contribution to the 
induced Newton constant and the cosmological constant. 

To obtain the suitable mass spectra, we use the method of dimensional deconstruction 
[2] and its generalization [3]. In the generalization of the deconstruction model based on 
a graph, the eigenvalues of the graph Laplacian and the adjacent matrix gives the mass 
spectrum of the particle. Thus we can easily control the induced quantities at one-loop level 
in such a model [4]. 

Wc also study self-consistent static solutions for a static Einstein universe in a graph-based 
induced gravity. We have considered self-consistent Einstein universe at finite temperature 
in [5] . In the present paper, we use the calculation method with the spectral density function 
of the graph and search for the static solution supported by the degenerate pressure of the 
fermion at zero temperature. 

The present paper is organized as follows. In §2, we will examine the UV-divergences 
in field theory with the heat kernel method. The way to construct suitable models using 
the knowledge of the graph structure is shown in §3. In §4, divergences in the effective 
gravitational action are regularized for a static Einstein space. It is shown that the technique 
with the density function to evaluate the effective action for an Einstein space in §5. In §6, 
strongly-degenerate fermions and a self-consistent solution in our model is studied. We give 
a summary and future prospects in the last section. 
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II. UV-FINITENESS CONDITION 

Induced gravity has been studied by many authors [1] . In terms of the heat kernel method 
[6], the one- loop effective action can systematically be expressed as an integral form using 
Schwinger's proper time. 

The classical action for a free field can be written as a quadratic form with a differential 
operator on the spacetime manifolds. The operator trace (Tr ) can be evaluated by the 
standard way to rewrite 

ilVln// = -i/^ fiv[e-^], (2.1) 

where H is a. Hessian operator which appears in the free-field action. The heat-kernel 

expansion can be expressed as, in four- dimensional spacetime, 



Tr 



^ J d^x^jdet^ [tr oo + nr oi + t'^ tr 03 + o{t^)] , (2.2) 



where Qfj^^, denotes the spacetime metric and tr means the trace over the spacetime indices. 
The Seeley-DeWitt coefficients Gp (p = 0, 1, 2, . . .) depend on the background fields and the 
first few coefficients have been known for several types of wave operators. The one-loop 
effective action for the background fields is given by the collection of the contribution of 
various matter fields to the heat-kernel coefficients. 

It is straightforward to see where the UV divergences occur, which we are interested in. 
The UV divergences arise from the integration in the vicinity of t = 0. These divergences 
arise from the first few terms of the heat-kernel expansion. If we manage to introduce a 
UV-cutoff scale A, the lower bound of the integration on t is replaced to 1/A^. To seek 
the condition for cancellation of UV divergences from various matter fields, we need only to 
consider massless fields. In the present paper, minimally-coupled scalar fields, spinor fields, 
and vector fields are taken into consideration. 

The first Seeley-DeWitt coefficient ag, which is a constant value, has been found for such 
fields. The value for each mode is: ao = 1 for a scalar mode, ao = 2 for a spinor field, and 
ao — 2 for a massless vector. Then the effective Lagrangian at one-loop level includes the 
following cutoff-dependent term proportional to [1] 

-L^{No-2Ny, + 2N,)A\ (2.3) 

where A^o is the number of minimal scalar degrees of freedom, N1/2 is the number of two- 
component fermion fields, and A^i is the number of massless vector fields. Note that the 
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spinor field contributes with a negative sign for its fermionic nature. The expression (2.3) 
corresponds to the cosmological constant or dark energy, if we treat it as a cutoff-regularized 
theory. 

The less divergent term comes from the coefficient ai. The coefficient for each mode is: 
Oi = R/Q for a scalar mode, Oi = —R/6 for a spinor field, ai = —2R/3 for a massless vector 
field, where R is the scalar curvature of the spacetime. Thus the coefficient oi leads to 
the induced Einstein- Hilbert term. The effective Lagrangian at one-loop level includes the 
following cutoff-dependent term proportional to [1] 

1 



1927r- 



:iNo + Ni/2-mWR- (2.4) 



Now we find that, to cancel the quartic and quadratic divergent terms, which diverge as 
A — > oo, we should choose 

No^2N, Ni/2^2N, Ni^N, (2.5) 

where = 1,2,3,.... 

The value of the Seeley-DeWitt coefficients can be confirmed when we set a specific 
background-space geometry. Since the eigenvalues of the wave operators for various fields 
on are well known, the trace part {tr cxp [—Ht] } for each field can be evaluated as follows 
and has an asymptotic form for small t: [5] 



oo 



a2 



5^2£(£+ l)exp 

oo 

5]2(£2-l)exp 
e=2 



(£+1/2)2^- 



a2 



(47rt)3/2 

2(27r2a3) / 1 



1 + -^t + ■■■], for a scalar mode , (2.6) 
J 



(47rt)3/2 V 2a2 
2{27r^a^) ( 2 



1 — — + • • • , for a spinor field , (2.7) 



(47ri)3/2 V a? 



1 -t + ■ ■ ■ , for a massless vector , (2.8) 



where a is the radius of -S"^. Because we know that the volume of is 27r^a^, the scalar 
curvature of is Q/a? and the second-order (Euclidean) time-derivative contribution gives 
a factor (47rt)~^/^, the values of oq and ai for these fields mentioned above can be verified 
by (2.6-2.8). 

If the massless matter content satisfies the condition (2.5), there is no quartic nor 
quadratic divergence and also no induced gravitational action because of absence of mass 
scales. Thus we should consider massses of the fields to yield the finite contribution of 
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quantum effects. Nonetheless, for cancellation of UV divergences, the condition (2.5) is still 
necessary. 

The algorithm to include the masses is very easy in the Schwinger time integration. We 
only attach the following to the integrand for each field 

J2 e-(-^)^* ^N,-t j:{m% + f- j:im% + ■ ■ ■ ^ - tTr + t' -Tr + ■ ■ ■ , C^-^) 

1=1 i=i ^ i=\ ^ 

where is the mass-squared matrix for spin-s field. 

In addition we need some interpretations in this trick. For massive spinor fields, we replace 
A^(y2 spinor fields to massive N'^j^l'^ Dirac fields. For massive vector fields, we replace N'^ 
massless vector fields as transverse modes and N'^ scalar modes as longitudinal modes to N'^ 
massive vector fields. Now we find the additional quadratic divergence is proportional to 

TrM| -4TrM|, + 3TrM^, (2.10) 

where M| is the mass-squared matrix of A'^^ massive scalar fields, M|, is that of iV(y2/2 
massive Dirac fields, and My is that of N'^ massive vector fields. 

Finally, the condition for cancellation of the quartic and quadratic divergences is con- 
cluded as follows. The matter content is: IN — Nq — N[ massless scalar fields, 2N — N[^2 
massless Weyl spinor fields, N — N[ massless vector fields, A^q massive scalar fields, N'l^^/'^ 
massive Dirac fields, and N[ massive vector fields. Moreover, massive fields must have mass 
matrices which satisfy Tr M| - 4 Tr M|, + 3 Tr = 0. 

III. GRAPH-BASED CONSTRUCTION OF A SPECIFIC MASS MATRIX 

In this section we construct the field theory with suitable mass matrices which satisfy the 
UV-finite condition expressed in the previous section. 

Now we remember the concept of dimensional deconstruction [2] , which is equivalent to 
considering a higher-dimensional theory with discretized extra dimensions at a low-energy 
scale. A moose diagram is used to describe this theory, and is no more than a graph. The 
A'-sided polygon is identified as an example of simple graphs, a cycle graph Cjv. 

A graph G consists of a vertex set V and an edge set £, where an edge is a pair of distinct 
vertices of G. The degree of a vertex v, denoted by deg{v), is the number of edges incident 
with V. If all the degrees of vertices of a graph are equal, we call such a graph as a regular 
graph. 
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We can consider the orientation of an edge. The graph with directed edges is dubbed as 
a directed graph. An oriented edge e = [u, v] connects the origin u — o(e) and the terminus 
V = t{e). 

Spectral graph theory is the mathematical study of a graph by investigating various 
properties on eigenvalues, and eigenvectors of matrices associated with it [7]. Now we 
introduce various matrices that are naturally associated with a graph [3, 7] for later use. 

The incidence matrix E{G) is defined as 



The adjacency matrix A{G) is defined as 

{■A)vv' — 

The degree matrix D{G) is defined as 



1 if f = o(e) 
-1 iiv = t{e) 
otherwise 



1 if V is adjacent to v' 
otherwise 



deg{v) ilv — v' 
otherwise 



(3.1) 



(3.2) 



(3.3) 



Note that Tr A = and Ti = TiD, and for a regular graph, D is proportional to the 
identity matrix. 

The graph Laplacian (or combinatorial Laplacian) A(G) is defined as 

deg{v) iiv — v' 

—1 if V is adjacent to v' ■ (3.4) 
otherwise 



(A)„„, = {D- A\ 
The most important observation is 



(3.5) 



A = EE^ , 

where E^ is the transposed matrix of E. The Laplacian matrix is symmetric, so its eigen- 
values are non-negative. Note also that Tr A = TrD and Tr A^ = TrD^ + TrL>. 

The simplest model of vector fields has been studied by Hill and Leibovich [8]. The 
generalized model associated with a general graph is written down as [3] 



(3.6) 



vev 



eeS 
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where the covariant derivative is 

P^C/e^(9'^ + iA^'(,)-i>l^(,))C/e, (3.7) 

with \Ue\ — f, f is a constant with the dimension of mass. The vector fields A!^ are assigned 
at vertices of G and the scalar fields Ue are assigned at edges of G in this model. 

Similarly, any kind of fields can be associated with a graph and their mass-squared matrix 
can be written using the graph Laplacian. For scalar fields, we assign a scalar field 0^ to 
each vertex v of G. A difference can be defined on each edge e as 

d4>e = 4>t{e) - (t^oie) = " X] ^ev(t>v ■ (3.8) 

Thus a mass term for scalar fields can be constructed as 

PY.d^edct>e = fY. E <Pv'E,>eEjj, = f ^ <Pv^vv'<Pv' ■ (3.9) 

ee£ ee£v,v'ev v,v'ev 

For spinor fields, the mass term can be expressed using the incidence matrix E. For 
example, the Lagrangian density of fermion fields can be written as [3] 

-Y^iljRypi/jR, -Y^^LeP^Le - / ^ E [(^ie(^'^)e.^R. + h.c] , (3.10) 
veV eeS ee£ veV 

where the subscripts L and R denote left-handed and right-handed fermions, respectively. 
Namely, the left-handed fermions are assigned to the edges while the right-handed ones are 
assigned to the vertices. The mass-squared matrix for ipny is expressed as pEE^ — /^A 
while that for ipLe is pE'^E = pA. The matrices A and A have the same spectrum up to 
zero modes. Thus the mass spectrum of fermions governed by the Lagrangian (3.10) is also 
given by the eigenvalues of the graph Laplacian (3.5). For details, see Ref. [3]. 

With the knowledge in spectral graph theory [7] , wc can find that the UV divergent terms 
are concerned with the graph Laplacian. Therefore, the UV divergences can be controlled 
by using the graph Laplacian and we can construct the models of UV-finite induced gravity 
from spectral graph theory. 

A prescription is as follows. First we prepare three graphs, Gs, Go and Gy- All these 
graphs have N vertices. We can construct Lagrangians whose mass-squared matrices satisfy 

TrM| = TrM^ = TrM^ , TrM| = TrM|, = TrM^ , (3.11) 
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^ Tr [e-^i* - 4e-^^* + Se"^^*! , (3.12) 



1 ^ Tr fe-^i* + 2e-^^* - Se"^^*! . (3.13) 

Ltt)^ Jo t'^ J 



by choosing graphs as D{Gs) — D{Gd) — D{Gv) [4]. Then we find that the induced vacuum 
energy at one-loop level is [4] 

1 /-"^ dt 

and the inverse of the Newton constant is given by [4] 
1 

IGttG ~ ~ 6(47 

In the flat-space limit, the one-loop vacuum energy has been calculated for fleld theory 
associated with the cycle graph C„ [4] . The degree matrix of a cycle graph C„ is an n x n 
diagonal matrix diag.{2, 2, ... ,2). We select a type of non-simply-connected graphs G{nj} = 
C„^ UC„2 U- • • — \Jj2 ni=N Cfii: which has N vertices. The degree matrix of G^m} is an NxN 
diagonal matrix diag.{2, 2, ... ,2). Therefore, if the mass-squared matrix is proportional 
to the graph Laplacian of TrM^ and Tr (M^)^ are independent of the choice of the 

set {rii}, as long as J^i^i = N is fixed. We can choose different sets {rii} for scalar, Dirac, 
and vector field model in order to obtain non-zero value for the Newton and cosmological 
constants [4]. 



IV. EVALUATION OF THE EFFECTIVE ACTION IN WITH ZETA 

FUNCTIONS 

We will consider a model for the static universe with spatial topology 5"^ with the radius 
a, in later sections. The self-consistent induced gravity model at finite temperature has 
been studied in Ref. [5]. We will study degenerate fermions at zero temperature and the 
self-consistent universe later in the present paper. 

In this section, we evaluate the one-loop vacuum energy for the spacetime R x S^'. To 
this end, we use (2.6-2.8) in the Schwinger integral form of the effective action. Here we first 
integrate over the proper-time t, but then we slightly shifted the power of t in the integrand. 
For example, an expression which appears in the effective action is rewritten as 



r(s-i/2)^ (i + l) 



til + 2) 2 ' 



a? 
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We then look for where divergences occur. We follow an analogous method used in Ref. [9], 
to separate a convergent summation from others. Now we convert it to 



Eofm^a^) = — I/O — / I/O 

t'o W + 2) + m2a2]^-^/2 + w?o? - l]'"'/' 

e 



E 



+a(2s - 3) + (1/2 - s)(m2a2 - l)C«(2s - 1) + ^ ^Cii(2s + 1) , (4.2) 



where C,e{z) is the Riemann's zeta function. Similarly we find 



£(£+!) 



E£)(m^a^) = 

[(^+1/2)2 + m2a2]^-^/' 1/2)2 + ^2q2]--V2 



4E 



(£+1/2)2-1/4 



-4E 

^=0 



(£ + 1/2)2-1/4 
[(£ + 1/2)2 + m2a2]*-'/' 
(£+1/2)2-1/4/^^ (l/2-g)m2a2 (s2 - l/4)mV^ 



(£ + 1/2)2^-1 



(£+1/2)2 



2(£+ 1/2)^ 



+4U2'^-'-1)Cr(2s-3) + 



2 2 

s I m a 



(2^-1 - 1)C^(25 - 1) 



+ 



s2-l/4)mV (l/2-s)m2a2 



(s2-l/4)mV 



(2^^+^ - l)C«(2s + 1) 



.(22^+^-l)C«(2s + 3) , 



(4.3) 



and also 



Ey(mV) = 2^^ 



£2 



2E 



ti [£2 + m2a2]^-'/2 ^tt + m2a2]^-'/' 

'2 - 1 £2_^ / (i/2-g)m2a2 (g2 _ i/4)mV ^ 



2E 

€=1 



[£2 + m2a2]^-V2 £2.-1 
1 



+2 Ci?(2s - 3) + 



+ 



l/4)mV 



s ) m^a^ — 1 

2 



2 2 

s\m a 



£2 

a(2s - 1) 
Cii(2s + 1)- 



2£4 



l/4)mV 



■Cii(2s + 3)| . 

(4.4) 



Since Cii(— 3) = ^io Cii(~l) = finite, only divergent part for s — >■ in each E is 

the term including Cii(2s + 1). The divergent parts are 



2Ef (m^a^) = (^2 _ 1^ ^^2^2 _ X)\^{2s + 1) , 



(4.5) 



Eg^(m'a') 



2s^ - - 1 m^a^ - ( - - s 1 wfa 



,2„2 



(2^^+^ - l)C«(2s + 1) , (4.6) 
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4 4 

m a — 



(1 - 2s) m^a^ Ci?(2s + 1) • 



(4.7) 



In the graph-based model reviewed in the previous section, we can set Tr M| = Tr = 
Tr My as well as Tr M| — TrM|, — Tr My. Thus the divergence in the induced action is 
proportional to 



This residual divergence is only in E5 and independent of mass, in other words, it appears 
even in the case with massless (minimal) scalar fields. Ehzalde [10] argued that this diver- 
gence should be dealt by 'principal part prescription'. In the prescription, the pole term in 
the Riemann's zeta function is discarded. This minimal subtraction yields 



where 7 is the Euler-Mascheroni constant (7 ~ 0.577216). 

Apart from the divergence, the divergent terms up to have been canceled. Corre- 
sponding to the analysis by using integral form ala Schwinger, we find that divergences 
including mass parameter can be cancelled in our graph-based models. 

We now redefine the finite part of summations as 




Cr(2s + 1) = ;^ + 7 + 0(s)^7, 



(4.9) 



E'5(mV)=5]£2 V£2 + m^a^ - 1 - £ 1 + 




) 



(4.10) 



00 r 

J:'^{m'a') = 4 5: [(£ + 1/2)2 _ y^]^ ^(£+1/2)^ + 



m?a? 



i=0 




(4.11) 



and 




(4.12) 



Then we find the effective action in the form. 



(4.13) 
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V. USE OF SPECTRAL DENSITY FUNCTION OF A GRAPH 

In this section, we introduce the spectral density function of a graph [11]. The use of the 
spectral density makes the analysis of the Casimir energy very easy. In the present paper, 
we consider only regular graphs. Remembering that the graph Laplacian is expressed as 
A = £) — A, we need only to consider the spectral density function for the adjacency matrix 
A in the case with a regular graph. 

We start with the case for a cycle graph Cn, for example. The spectrum of the eigenvalues 
for the adjacency matrix of is 

077- i- 

Afe = 2cos^, (A; = 0,1,...,7V-1) (5.1) 

and thus the eigenvalues for A are = 2 — 2 cos ^ = 4sin^ It has been shown [11] 
that, since 

li^ ^ E /(A.) = f'f{2cosnt)dt = - f f{x)^^=, (5.2) 
the spectral density in the large N limit can be employed as 

/+00 \ r2 (^2; 

f{x)pN{x)dx = - f{x)-j==. (5.3) 
-00 TT J-2 V4 — X^ 

Namely, the summation about the discrete eigenvalues becomes an integration over the 
continuous variable x with the spectral density function Poo{x), 



1 1 



for - 2 < X < 2 



TT V4 x'^ fQj. cycle graphs , (5.4) 

otherwise 



in the large N limit. Incidentally, the precise spectral density function for Cat with a finite 
N is known as 

I otherwise 

where Tn{z) denotes the Chebyshev polynomial. 

The spectral density function is known for other several graphs. The trace formula for 
regular graph G of degree g + 1 on A'" vertices is [12] 

- f e*^' - /'^ e-^^i^dx + 1 V V -^h. .{2 rat) (5 6) 

~ 27r J-2^,^ (5 + l)2-x2^'^+ Ar4-;^2W2^^^(^)^V?^), (5.bj 



g k=i 
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where g runs over the set of all oriented primitive closed geodesies in G, and i{g) is the 
length of g, while In{z) is the modified Bessel function of the first kind. Then 

q+l \/4g- 



Pooix) 



^-nnw^ for - 2^<x <2^ . . x , ^ 

^ V V for + 1)— regular graphs . (5.7) 

otherwise 



In the present paper, we will concentrate ourselves on the case with the graph G^m} = 
Gni U C„2 U • • • = [jj2.ni=N^ni- Clearly enough, one find that the spectral density function 
Poo is independent of the choice of {rij}. 

This fact imphes that the finite contributions for the Newton and cosmological constant 
come from the pn — Poo if the summation is evaluated as the integration over the continuum 
variables. Therefore the Casimir energy behaves as 1/a^ x 27r^a^ and the similar contribution 
which dominates if a is small are substantially calculated only by using poo and that is 
independent of values for the Newton and cosmological constant in the flat-space limit. 
This universal conclusion may be interesting if we try to extend the present approach to the 
case with general graphs. 

Turning to the present analysis, we assume that the mass-squared matrix is given by 
pA{G), where / is a unique mass scale in the model. For large iV, the effective action, 
where the Casimir energy is dominant, becomes 



ma) ^ Ya 1-2 [^'^^^'"'(2 - x)) - T!j,{fa\2 - x)) + i:'y{fa\2 - x)) 



N 



ttVT— X 



=dx. 
(5.J 



In the next section, using this result, we study a self-consistent cosmological solution for 
an Einstein universe in the graph-based induced gravity model. 



VI. DEGENERATE FERMIONS AND A SELF-CONSISTENT UNIVERSE 

We consider a model for the static universe with spatial topology with the radius 
a. The static homogeneous, closed space is often called an Einstein universe. The self- 
consistent induced gravity model at finite temperature has been studied in Ref. [5]. In the 
present paper, we study the self-consistent cold universe at zero temperature and we will 
consider degenerate fermions. Although the cold universe seems to have less relevance to the 
actual universe than the hot case, it can be a possible phase between quantum cosmology 
and classical cosmology. 



13 



In the static spacetime, it is known that the effective action can be interpreted as the 
total free energy of the quantum fields at finite temperature [13] . Similarly, we consider the 
thermodynamic potential for the case with a finite chemical potential. 

The thermodynamic potential of a system of strongly-degenerate fermionic fields at zero 
temperature can be computed as [14] 



127r2 



5 2 
2^. 



/i 



+ ^rrii In h 



, (6.1) 



where /j, is the chemical potential and 9{y) is the step function, 9{y) ~ 1 for y > and 

ei^y) = for y < 0. 

For the case with the model associated with the graph which consists of a set of Cn, 
in the large N (the total number of vertices) limit can be reduced to 



X 



X 



127r2 



6{n — m{x)) 



li\J — m?{x) — 
dx , 



-m'^fx) In 
2 ^ ' 



/i 



m[x 



) \ 



m^ix) 



N 



(6.2) 



with m'^{x) = p{2 — x). 

It is known that the fastest way to obtain self-consistent equations is by using the total free 
energy in the finite-temperature case [15]. Similarly, we consider the total thermodynamic 
potential Q as the sum of the contribution of quantum effects Qq derived in the previous 
section and that of degenerate Dirac fields Vlu . The energy of the system is given by 



where 



(6.3) 



(6.4) 



is the fermion number, which suffers no correction from Qq. The pressure P is obtained by 



Px(2.V) = -3a^, 

as in the finite-temperature case. 

The self-consistent equations can be derived as 



(6.5) 



(6.6) 
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and 

where the first equation corresponds to the 00-component of the Einstein equation with one- 
loop corrections and the second corresponds to the diagonal component in a spatial direction. 
Thus the extremal point of ii~^fl{fa, f / ii) provides a solution to the self-consistent equation. 

1 

. 8 
0.6 
. 4 
0.2 



0.2 0.4 0.6 0.8 1 

a 

FIG. 1. A contour plot of ^n"^^. A solution of the self-consistent equation can be found at the 
maximum point. 

In FIG. 1, we show the contour plots for fi/// obtained by numerical calculations, whose 
extremum provides a self-consistent solution. The horizontal axis indicates the scale factor 
a, while the vertical one l//i, in the unit of /. 

Since the Casimir energy is dominant for small a, the solution can be found at the 
maximum of corresponding to the Casimir regime defined in Ref. [15]. The stability 

is not expected, for the extremum of the potential is actually the maximum point. 

VII. SUMMARY AND OUTLOOK 

In the present paper, we have examined ultra-violet divergences of a one-loop calculable 
model for induced gravity. We have found that finite values for the Newton and cosmological 
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constant can be realized if the mass-squared matrices for scalar, spinor, and vector fields 
satisfy a few conditions. 

It has been found that the model which has the suitable mass matrices can be obtained 
by the graph-based construction. In this paper, we focused on a type of the regular graph 
such as G = Cm U U ■ ■ ■. 

To evaluate the effective action for an Einstein universe, we need the knowledge of graph 
spectrum. We have introduced the spectral density function of the graph and found that it 
is useful to calculate the Casimir-energy dominant case, for small a and large N. 

The spectral density is also convenient to evaluate the thermodynamical potential of 
strongly-degenerate fermions. We have studied self-consistent Einstein universe at zero 
temperature with degenerate fermions in our model. We found that the Casimir regime can 
been seen. 

In the present analysis, we have constructed models using cycle graphs, but we are also 
interested in the model of general graphs. As future works, trace formula for a regular graph 
[12] will be useful. 

The universal behavior of the effective action for large N and small a under the condition 
of the fixed type of the associated graph, is interesting. If the construction of the model with 
dynamical selection of graphs is possible, say, utilizing the Higgs-like mechanism assigned at 
edges or vertices, it can be imagined that many large-scale universe with different Newton 
and cosmological constants would develop once from a single state with a large Casimir 
energy. Anyway, we should investigate some variation of the present model. 
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